7 

AD-A119  997  PUERTO  RlCO  UNIV  RIO  PIEDRAS  DEPT  OF  PHYSICS  F/6  7/4 

THE  APPLICATION  OF  THE  MODIFIED  60UY-CHAPMAN  THEORY  TO  AN  ELECT— ETC (U) 
•  SEP  82  L  B  BHUIYAN'  L  BLUM*  D  HENDERSON  N00014-81-C-0776 

UNCLASSIFIED  NL 


AD  A119997 


t 

i 


OFFICE  OF  NAVAL  RESEARCH 
Contract  N001-81-C-0776 


Task  No.  NR  051-775 


The  Application  of  the  Modified  Gouy-Chapman  Theory 
to  an  Electrical  Double  Layer  Containing  Asymmetric  Ions 


by 

L.  B.  Bhuiyan*,  L.  Blum*,  and  D.  Henderson** 


Prepared  for  Publication  in 
The  Journal  of  Chemical  Physics 


*  University  of  Puerto  Rico 
Department  of  Physics 
Rio  Piedras,  Puerto  Rico 


**  IBM  Research  Laboratory 
San  Jose,  California 


September  21,  1982 


Reduction  in  whole  or  in  part  is  permitted  for 
•i  any  purpose  of  the  United  States  Government. 


*This  document  has  been  approved  for  public  release 
and  sale,  its  distribution  is  unlimited. 

Q_ 

CZ> 

< 

*This  statement  should  also  appear  in  Item  10  of  Document  Control  Data  -  DD 
_ I  Form  1483.  Copies  of  form  available  from  cognizant  contract  administrator. 


Copy  available  to  DTIC  does  not 
permit  fully  legible  repioduction. 

82  10  07  024 


I 


✓ 


<CCUfttTV  CLASSIFICATION  Of  THU  PAGE  Onto  Intorop 


REPORT  DOCUMEHTATION  PACE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

i.  REPORT  NUMBER  j.  OOVT  ACCESSION  NO. 

Technical  Report  #  '  3  &£>  -  A//  9  90  7 

i.  recirient-s  catalog  number 

A-  TITLE  (on*  SuBtiUo) 

The  application. of  the  modified  Gouy-Chapman  Theory 
St  a  an  eieStfical  double  'layer  containing  asymmetric 
"ions."  *w  ■  - 

S.  TYRE  or  REPORT  A  PERI  00  COVERED 

Interim  Technical  Report 

S.  PERFORMING  orc.  report  NUMBER 

7.  author**;  _  . 

L;  B.'Bhuiyan,  L.  Blum  and  D .-  Henderson 

s.  contract  or  grant  number**; 

N0001 4-81-C -.0.7  76 

*•  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Department  of  Physics 

University  of  Puerto  Rico 

Box  AT,  Rio  Riedras,  PR  OCmi"  *  .*  ' 

It.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS  ^ 

Code  472 

Office  of  Naval  Research 

Arlington,  Va.  22217 

ta.  report  date 

June,  1982 

13.  number  of  pages 

18 

I*.  MONITORING  AGENCY  NAME  A  AOO H£SS(tt  alttormni  from  Controlling  OUlcm) 

IS.  SECURITY  CLASS.  *•*  Olio  roport) 

unclassified 

OISTmTSutTon  STATEMENT  (ml  u Us  TtZpmrt) 


Approved  for  public  release;  Distribution  Unlimited 


17.  DISTRIBUTION  STATEMENT  (ot  Him  atNTMl  mntmrmH  In  Block  20.  II  dlllmrmnt  tea  R.port} 


>».  SUI*»LSMtHT  ARY  NOTES 

Prepared  for  publication  in  Journal  of  Chemical  Physics 


19.  KEY  WORDS  (Continue  «n  r«v»n«  ud*  1/  n«ct«««y  <nd  Identify  by  block  numbmr) 

Gouy-Chapman ,  electrode,  interface,  asymmetric  ions. 


vO.  aCSTSaCT  'Continv  on  to+mt*o  #*tf»  l i  /*•€»•* ory  md  leontify  by  block  ‘weiptr; 

^"The  modified  Gouy-Chapman  Theory  is  applied  to  1:1,  2:1  and  1:2  electrolytes 
of  unequal  sizes  near  a  flat  electrode.  This  mode^  features  a  potential 
of  zero  charge  in  the  absence  of  specific  adsorption. 


DO  I  JAN *71  1473  ESIT1CN  or  I  NCV  *3  IS  obsolete 
S/N  6»C2.L?-CW-d4Cl 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  BEST  QUALITY 
PRACTICABLE.  THE  COPY  FURNISHED 
TO  DTIC  CONTAINED  A  SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


INTRODUCTION 

The  accumulation  of  ions  near  a  charged  electrode,  resulting  in  a  so-called  electrical 
double  layer,  has  been  the  subject  of  considerable  recent  interest.  Both  integral  equation 
methods1*5  and  Monte  Carlo  simulations6  have  been  applied  to  this  problem. 


The  classical  theory  of  the  electric  double  layer  is  the  Poisson-Boltzmann  theory  of 

'  '  '  "  '  f  ‘  ‘  '  ■ 

'  '  '^  j'Gouv'  and  Chapman6  in  which  the  solvent  is  assumed  to  be  a  uniform  dielectric  medium  _ 

dielectric  constant  "is  e.  -  iv-  ***•  *  -  &  f 
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.  .  In  the  Gouy-Chapman  (GC)  theory,  the  ionxMesare  neglected.  It  is  a  theory  for 

-• */-*■.  •  •  -♦  • 

point  ions.  However,  if  the  ionic  profiles  are  displaced  by  a/2,  where  a  is  the  ionic  ^  ' 

.«* .» V..  .  _  ;■  .  „ .  . 
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diameter,  so  that  the  distance  of  closest  approach  is  a/2  rather  than  zero,  the  GC  profiles 


for  monovalent  ions  are  reasonably  good.  We  refer  to  this  modification  as  the  modified 

■  - -  . 
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J  Gouy-Chapman  (MGC)  theory.  In  the  MGC  theory,  the  hard-core  interactions  among  the 
r  "'i/L' v  ions  are  neglected,  but  the  hard-core  interaction  of  an  ion  with  the  electrode  is  not 

■VSr^- •;  .  -  — — : 

;'.>^,S^!^^^.neiglecte’dr'.The"MGC  is  less  satisfactory  for  multivalent  ions  but  is  still  of  interest  because 
-  it  is  a  useful  starting  point  for  more  accurate  interactive  procedures  based  on  integral  ,  ... 

;  -  v  .  .  .  ......  .  '  . 

-  equations.  .  V  .  -  -  ••*••* 

. :  The  GC  theory  was  originally  applied  to  symmetric  ions.7,8  'ine  extension  to  2-1  and 

1-2  electrolytes  was  obtained  by  Grahame.9  Recently,  Valleau  and  Torrie 10  have  applied 
the  MGC  theory  ions  which  are  symmetric  in  charge  and  asymmetric  in  diameter.  In  this 
note,  we  apply  the  MGC  theory  to  ions  which  are  asymmetric  in  both  diameter  and  charge. 
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THEORY 

/  The  starting  point  of  the  GC  is  the  one-dimensional  Poison-Boltzmann  equation 


d2* 


—■  -  2  ziPi  exP  {0zie<P(x)} 

dxz  e  i-l 


0) 


where  <p(x)  is  the  mean  electrostatic  potential  at  a  distance  x  from  the  electrode,  e  is  the 


temperature.  The  parameter  m  is  the  number  of  species  of  ions  in  the  region  of  space  to 

l  .v  •  .  *C:.  •.  . 

^.^'-^-^.;  :t-''which  Eq.  (1)  is  being  applied.  If  m»0,  the  RHS  of  Eq.  (1)  is  zero  and  Eq.  (1)  becomes 
Laplace’s  equation  for  charge  free  space,  “*&.***:  ./  .  r.  v.: 


>  :  Independently  of  the  value  of  m>  1,  Eq.  (1)  can  be  integrated  to  give 
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T‘  where  A  is  a  constant  of  integration.  ....  "  ^  .L 

Ja»S5(WtTytS**^’v,r.*  w '  •  ' '  '  ■/'  '  ’  •  '-■'  ■  '■  '  1 t*?1 
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-  'i  ^  _  To  obtain  $(x),  Eq.  (2)  must  be  integrated.  This  can  always  be  done  numerically. 

t..*.  i-  t  However,  results  in  closed  form  have  been  obtained  only  for  m->0,  1  and  2  with  either  all 

vvv*'"  '■  •••■■■ 

the  Zj  equal  (symmetric  charge)  or  with  Zj-1  and  Z2«2,  or  equivalently,  Zj—2  and  z2-«l. 
Therefore,  we  restrict  ourselves  to  n-n,  1-2  and  2-1  binary  solutions. 


The  geometry  of  a  general  binary  system  is  shown  in  Fig.  1.  There  are  three  regions 
to  consider:  0<x<Oj/2,  and  x-o2/2,  where  <*j<a2  am*  a2  are  ^*e 

diameters  of  the  two  species  of  ions.  In  these  three  regions,  there  are,  respectively,  no  ions 
present,  only  ions  of  one  species,  or  both  species  of  tons. 


,.l  In  the  first  region  (m«0),  0sx^O|/2  and  integrating  Eq.  (2)  we  get 

*00  -  *<o)  +  x*'(0) . 

, .  v v;>;  where  *'(0)  is  the  derivative  of  *(x)  at  the  electrode.  In  the  second  region  (m-1), 

>  \  0\/2<tx<to2/2, 

Q  -Tv?  exo  fflz,e6fxU  +  a1*/2  . 


^^^^^gwhere  yj  «8irkT p [/ 1. ■  Integrating  gives  is.  $«$*-****v "w  ■*?!?* 

T  *:v it  -*  •  ..  .'•  K'-s~ 

/"/-I  a  T/  i+b  exp{-p2le*(x)}  x2  ,11  ... 

-m  B^L(r-B.xP{-^(x,))  -*Jj  •  A>0 

:  fcit?(x)  -<  ..  ..  )  -  (5) 

r?:?r.^'-ir?v':l  ^n{  +  tan  2(B_ x^Zjev^A/2)]!  ,  A<0  , 
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where  B  is  a  constant  of  integration. 
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golution  jor*(x)  in. the  third  region  (m-2),  x><72/2.  is  well  known.  For 
^^^^flrmmetnc^ ^charge,' 1  Zir«zl^ ^the  result  is7>8 
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For  the  unsynunetric  2-1  and  1-2  case,  the  result  is9 
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exp  [/5e*(x)]  - 
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with 
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1/2 

{2  exp  [/3e$(o2/2)]  +  1} 

.  ->  j *  ■  V; w  '  . 

*r  ^ .  4# '  1  /2 

{2  exp  [— j8e$((r2/2)]  +  1}  •„ 


In  Eqs.  (6)  and  (7),  y-x-flr2/2,  and 
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(10) 


In  obtaining  Eqs.  (6)  and  (7),  we  have  employed  the  condition  <£(x),  <f/(x)-*-0  as  x-*». 
.Equations  (7a)  and  (9a)  apply  to  the  case  where  the  electrode  has  a  positive  charge  and 


~  J^l8'  (7b)  and  (9b)  apply  to  the  case  where  the  electrode  has  a  negative  rharg* 


The  constants  of  integration  are  determined  by  requiring  that  $(x)  and  $'(x)  be 


continuous  at  x*Oj/2  and  oj/2  and  ^(0)  be  equal  to  the  applied  potential.  In  principle, 
analytical  expressions  could  be  developed  for  these  constants  of  Integration.  However,  in 
practice,  we  found  it  convenient  to  satisfy  the  boundary  conditions  by  an  iterative 
procedure. 


RESULTS 

The  potential  drop  across  the  double  layer  (measured  from  x*0)  is  plotted  in 
Figs.  2-5  as  a  function  of  the  charge  density  on  the  electrode.  As  would  be  expected  at 
large  charge  densities,  the  potential  drop  is  largely  determined  by  the  counterions  since  the 
co-ions  are  excluded  from  the  vicinity  of  the  electrode.  Consequently,  the  results  for 


positive  charge  densities  are  insensitive  to  the  value  of  ox/a2  since,  for  positive  charge 
densities,  the  counterion  always  has  the  diameter  ff2“4-2sA. 

•  For  negative  charge  densities  on  the  electrode,  there  is  a  wider  variation  in  the 

•  -** *».'  v  -  "  ' 

.  _  •  potential  since  the  counterions  now  have  a  variable  diameter  a1  <a2.  At  large  negative 

JV  '■•Sj,'** .  V.  ..  . . 

char*e  ^osities,  the  potential,  as  expected,  approaches  that  of  a  symmetric  electrolyte  of 
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r.  i  -  •  *T»*  .  '•  r- 


<■ 


.  .  ,  At  small  charge  densities  on  the  electrode,  the  smaller  ions  can  approach  the  wall 

more  closely.  This  results  in  a  separation  of  charge  for  ions  of  unequal  diameter  even  at 
zero  charge  on  the  electrode  and  a  resulting  nonzero  potential  difference  acrossithe  double 


yr .  -  -i.  . 


---  layer .  This  has  already  been  noted  by  Valleau  and  Torrie^®  for  ions  which  have  a 
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symmetric  charge. 


‘  -  The  charge  profile  and  concentration  profiles  for  a  O.SM  solution  with  <Tj/ct2—0.v  are 

.-v«i^crf'3vsMr' »»<<- >  -- -  '  ■  ■■■•  v  . . 

plotted  in  Figs.  6  and  7  for  the  case  where  the  potential  drop  across  the  double  layer  is  • 

‘ 

■'  -f  -■i~'  ••  . '  zero-  The  potential  profile  is  not  monotonic,  as  is  the  case  when  a^ma2  or  when  the 

i?'  •  v 

"  ;  .  r  "  potential  drop  is  appreciably  different  from  zero.  The  concentrations  profile  of  the  smaller 

.A  -  .  '  1  • 

— ;  /  ;•  ions  is  also  nonmonotonic  at  V*0.  ~  • 

As  is  seen  from  Figs.  8  and  9,  even  for  the  relatively  small  potential,  I V  |  *0.026 
volts,  the  concentration  profiles  are  rather  similar  to  those  for  the  equal  diameter  case.  At 
large  potential  drops,  the  concentration  profiles  near  the  electrode  will  approach  those  of 
the  symmetric  electrolyte  whose  charge  and  diameter  are  equal  to  those  of  the  counterions. 
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CONCLUSIONS 

The  MGC  theory  has  been  applied  to  electrolyte  which  are  asymmetric  in  charge  and 
size.  Within  the  MGC  theory,  the  most  interesting  results  are  obtained  in  the  region  of 
very  small  potential  drop  and  very  small  charge  density  on  the  electrode  because  at  large 
potentials  or  charge  densities,  the  co-ions  are  excluded  from  the  vicinity  of  the  electrode. 

.  Consequently,  the  counterions  dominate  and  the  double  layer  properties  approach  these  of 
v  a  symmetric  electrolyte  whose  charge  and  diameter  are  equal  to  those  of  the  counterion. 

If  the  ions  differ  in  size,  then  the  smaller  ions  can  penetrate  closer  to  the  electrode 

than  the  larger  ions.  At  large  charge  densities,  this  effect  would  be  negligible  compared  to  -  _  r 

) 

the  coulombic  effects.  However,  this  effect  is  significant  at  small  charge  densities.  In 
.  particular,  charge  separation  and  a  nonzero  potential  can  occur  even  when  there  is  not 
charge  on  the  electrode. 

It  is  conventional  to  interpret  experimental  results  in  terms  of  the  GC  theory  in  which 
_  ion  diameters  do  not  appear.  As  a  result,  noncoulombic  effects^uch  as  specific  adsorption,  ~ 
are  required  to  explain  such  phenomena  as  a  nonzero  potential  at  zero  charge  on  the 
electrode.  It  is  not  necessary  to  invoke  such  effects  since  unequal  ion  sizes  can  also  give 
rise  to  such  phenomena.  This  is  not  to  say  that  specific  adsorption  is  not  an  important  — 
phenomena  in  double  layer  studies.  However,  estimates  of  the  importance  of  specific 

adsorption  obtained  by  subtracting  equal  diameter  Gouy-Chapman  results  from 

_  *  •  .  -  - 

experimental  studies  may  well  be  misleading.  IVxeX'fe  CW'e.  VtyoJoW  cit 
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FIGURE  CAPTIONS 


Figure  1.  Double  layer  geometry. 

Figure  2.  Potential  drop  across  a  double  layer  as  a  function  of  the  charge  density  on  the 
electrode  for  the  case  where  Zj-1,  z2— I,  e-78.5,  T-298K,  a2» 4.25A  and  Oj/o  ~*1,  0.5 
and  0.1.  The  solid  and  broken  curves  give  the  MGC  results  for  concentrations  of  Vi  and 
0.05M,  respectively. 

) 

Figure  3.  Potential  drop  across  a  double  layer  as  a  function  of  the  charge  density  on  the 
electrode  for  the  case  where  Zj—1  and  z2=«-X  The  other  parameters  have  the  same  values 
as  in  Fig.  2  and  the  curves  have  the  same  meaning  as  in  this  figure. 

"  Figure  4.  Potential  drop  across  a  double  layer  as  a  function  of  the  charge  density  on  the _ 

electrode  for  the  case  where  Zj*>2  and  z2=-l.  The  other  parameters  have  the  same  values  as 
in  Fig.  2  and  the  curves  have  the  same  meaning  as  in  this  figure. 

Figure  5.  Potential  drop  across  a  double  layer  as  a  function  of  the  charge  density  on  the 
electrode  for  the  case  where  Zj«2  and  z2*-2.  The  other  parameters  have  the  same  values  as 
in  Fig.  2  and  the  curves  have  the  same  meaning  as  in  this  figure. 

Figure  6.  Potential  profiles  for  the  case  where  V»0  and  where  o2/o2«0.1  and  Zj»2,  z2«-l 
and  Z|*l,  z2«-2.  The  concentration  is  0.5M.  The  other  parameters  have  the  same  values 
as  in  Fig.  2. 


Figure  7.  Concentration  profiles  for  the  case  where  V-0  and  where  ax/a2  -0.1  and  Zj-2, 

Zj-O.l  and  Zj-1,  z2-2.  The  concentration  is  0.5M.  The  other  parameters  have  the  same 
values  as  in  Fig.  2. 


Figure  8.  Concentration  profiles  for  the  case  where  V-0.026  volts  and  where  aj/^-o.l 

Zl“ 1#  Z2”"2'  Theconcentration  is  0.05M;  The  other  parameters  * 
^ :'thave  the  same  values  as  in  Fig.  2.  -  •'!“  / 


. FlgUrC  9‘  Concentration  profiles  for  the  case  where  V— 0.026  volts  and  where  ol/o2*.Q.5 

-  and  Zj  — 1,  Zj— 2  and  z^— 2,  z^— -1.  The  concentration  is  0.5M.  The  other  parameters  have 
< \\-f. the  same  values  as  in  Fig.  2. 
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